ADAPTING THE TIME STEP TO RECOVER THE 
ASYMPTOTIC BEHAVIOR IN A BLOW-UP PROBLEM 
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Abstract. The equation ut = Au + uP with liomegeneous Dirichlet bound- 
ary conditions has solutions with blow-up if p > 1. An adaptive time-step 
procedure is given to reproduce the asymptotic behvior of the solutions in the 
numerical approximations. We prove that the numerical method reproduces 
the blow-up cases, the blow-up rate and the blow-up time. We also localize 
the numerical blow-up set. 



1. Introduction. 

We study the behavior of an adaptive time step procedure for the foUowing 
parabohc problem 

( ut^ Au + uP in X [0,r), 

(1) <^ u{x,t) = ondnx [0,T), 

[ u{x,0) — uo{x) > onQ. 

Where p is superhnear (p > 1) in order to have solutions with blow-up. We assume 
uq is regular and C M'' is a bounded smooth domain in order to guarantee that 
w e C^'^. A remarkable fact in this problem is that the solution may develop singu- 
larities in finite time, no matter how smooth uq is. For many differential equations 
or systems the solutions can become unbounded in finite time (a phenomena that 
is known as blow-up). Typical examples where this happens are problems involving 
reaction terms in the equation like where a reaction term of power type is present 
and so this blow up phenomenum occurs in the sense that there exists a finite time 
T such that limj^T OIloo — +oo for initial data large enough (see [SI m and 
the references therein). The blow-up set, which is defined as the set composed of 
points X G such that u{x, t) +oo &s t ^ T , is localized in small portions of 
n, in [51] is proved that the [d — 1) dimensional Hausdorff measure of the blow-up 
set is finite. The blow-up rate at these points is given by u{x,t) ~ (T — t)~'p^ , 
moreover 



(see II1,[IS1). 
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We remark that these resuhs hold if p is subcritical {p < if d > 3). For 
supercritical p the solutions may present different behaviors. For that reason we 
assume p is subcritical along the paper, however the asymptotic properties of the 
numerical schemes described above hold for every p > 1. This is a difference 
between the continuous solutions and their approximations. 

Since the solution u develops a singularity in finite time, we investigate what 
can be said about the asymptotic behavior (close to the blow-up time) of numerical 
approximations for this kind of problems. In 16 the authors analyze a semidiscrete 
scheme (keeping t continuous) in an interval, = (0, 1). Here we generalize that 
analysis for several space dimensions an introduce the adaptive discretization in 
time. 

For previous work on numerical approximations of 1^ we refer to HI, H], [7|) 
[S], jni) the survey ^ and references therein. 

As a first step to introduce the method we propose a method of fines, that is: 
we discretize the space variable, keeping t continuous. In this stage we consider 
a general method with adequate assumptions. More precisely, we assume that for 
every h > small {h is the parameter of the method), there exists a set of nodes 
{xi, . . . , xn} C (N — N(h)), such that the numerical approximation of u at the 
nodes Xfe, is given by 

U{t) = (ui(t),...,u^v(t)). 

That is Uk{t) stands for an approximation of u{xk,t). We assume that U is the 
solution of the following ODE 

(2) MU'{t) = -AU{t) + MUP, 

with initial data given by Uk{Q) — uo(xk)- In ||2Jl and hereafter, all operations 
between vectors are understood coordinate by coordinate. 

The precise assumptions on the matrices involved in the method are: 
(PI) M is a diagonal matrix with positive entries ruk- 

(P2) A is a nonnegative symmetric matrix, with nonpositive coefficients off the 
diagonal (i.e. a.y < if i 7^ j) and an > 0. 

(P3) Ef=i a.fc > 0. 

As an example, we can consider a linear finite element approximation of prob- 
lem Q on a regular acute triangulation of (see [S]). In this case, let Vh be the 
subspace of HQ{n) consisting of piecewise linear functions on the triangulation. 
We impose that the finite element approximation : [0, T/i) — > verifies for each 

t e [0, Th) 

I {{uh)tvY = - f VuhVv+ f {{uhrvY, 
Jn Jn Jn 

for every v € Vh- Here {■)^ stands for the linear Lagrange interpolate at the nodes 
of the mesh. These conditions imply that the vector U{t), the values of Uh{-,t) at 
the nodes Xk, must verify a system of the form (0). In this case M is the lumped 
mass matrix and A is the stiffness matrix. The assumptions on the matrices M 
and A hold as we are considering an acute regular mesh. We observe that in this 
case Uh — . 
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As another example, if is a cube, fl = (0, 1)"^, we can use a semidiscrete finite 
differences method to approximate the solution u{x,t) obtaining an ODE system 
of the form (j^Jl . 

We also need some kind of convergence result for the scheme, we will state the 
precise hypotheses concerning convergence in the statement of each theorem. Fi- 
nally, in the Appendix we prove an L°° convergence theorem under the consistency 
assumption an give some examples. Now we state the two convergence hypotheses 
that we may need. 

(HI) For every r > ||u - Uh||L~(Ox[o,T-r]) ^ as /i ^ 
(H2) \\u — Uh||ifi(si)(i) ^ as /i ^ for a.e. t 

Writing these equations explicitly we obtain the following ODE system, 

N 

(3) mfeu^(i) = -^afciUj(i) + mfeU^(i), 1 < fc < iV, 

1=1 

with initial data Uk{0) = UQ{xk)- 

In the authors studied these kind of schemes in one space dimension, ob- 
taining results similar to those stated below. 

Once obtained the ODE system, the next step is to discretize the time variable 
t. In the authors suggest an adaptive time step procedure to deal with the heat 
equation with a nonlinear flux boundary condition. They analyze explicit Euler and 
Runge-Kutta methods, however all these methods have to deal with restrictions in 
the time-step. In this work we first analyze an explicit Euler method and next we 
introduce an implicit scheme in order to avoid the time-step restrictions. We use 
= . . . for the values of the numerical approximation at time tj, and 
Tj — tj+i — tj. When we consider the explicit scheme, is the solution of 



(4) 

or equivalently, if we denote duj^^ = ■^{u-'^^ 



N 

m,du^^ = - ^ ak^ul + m,(u^)P, \ <i<N 

k—1 

u°^uq{x^), l<i<N + l. 

While for the implicit scheme is the solution of 

MW+^ = MW - Tj {AW+^ + M{W)P) 
^' U{0) = ul 

Note that the scheme is not totally implicit since the nonlinear source is 
evaluated at time P while the discrete laplacian {A) is evaluated a time P^"^ . This 
mixture makes the scheme free of time-step restrictions while the explicit evaluation 
of {U^ Y avoids the problem of solving a nonlinear system in each step. 
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Now we choose the time steps tj = tj-^-i — tj in such a way that the asymptotic 
behavior of the discrete problem is similar to the continuous one. We will fix A 
small and take 

_ A 

where 

N 

— ^ nikul 

k=l 

This choice for the time step is inspired by PP, where the authors use an adaptive 
procedure similar to this. They adapt the time step in a similar way but using the 
maximum (L°°-norm) instead of (L^-norm). In their problem the maximum 
is fixed at the right boundary node (i.e. ||f7-'||cx; — ''^tv+i)- In this problem, the 
maximum (the node k such that u], — ||?7-'||oo) can move from one node to another 
as t goes forward. So it is better to deal with since, for example, we can compute 
its derivative. Anyway, as we use a fixed mesh for the space discretization we can 
compare both norms. The motivation of this choice for the time-step is that, as 
will be shown, the behavior of is given by 

Hence, with our selection of tj we can obtain 

+ 1 ^ + Tj{w^Y +\^W^ + {j + 1)A, 

and we obtain the asymptotic behavior of , which is, as we will see, similar to 
the one for the continuous solution. 

When we deal with the explicit scheme we will also require 



(7) A< min —(w°)p. 

l<i<N an 

Then we study the asymptotic properties of the numerical schemes. We will say 
that a solution of Q (or ©) blows up if 



lim ||[/-' ||oo = oo, and T/j.a := < oo, 

7 — ^oo ^ — ^ 

we call Th^x the blow-up time. To describe when the blow-up phenomena occurs in 
the discrete problem we introduce the following functional : M. 

<i>h{U) = {AU,U) - {-^MUP+\ME), 
p+l 

where iJ = (l,l,...,l). The functional $;i is a discrete version of 

This functional characterize the solutions with blow-up in the continuous prob- 
lem: in JI], it is proved that u blows up at time T if and only if $(u)(t) ^ — oo 
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as i — > r. We prove a similar result for the discrete functional and this allows 
as to prove that if the continuous solution has finite time blow-up its numerical 
approximation also blows up when the parameters of the method are small enough. 

Next we study the asymptotic behavior for the numerical approximations of the 
solutions with blow-up and we find that they behave very similar to the continuous 
ones. In fact we find that if Uh^\ is a numerical solution with blow-up at time Th,x 
its behavior is given by 

max uj - (Thx-t^y^^'-P'^l 

l<i<N ' ^ ' 

We use the notation f{j) ^ g{j) to mean that there exist constants c, C > 
independent of j (but they may depend on h) such that 



cgU) < fU) < Cg{j) 
Moreover, the numerical schemes recover the constant Cp in the sense that 

lim max ui{Th,x - i^y/(p-i) = Cp. 

j^oc l<i<N * ' ' ^ 

The functional $/i is also useful to prove convergence of the blow-up times. 
Unfortunately we can only prove the convergence of an iterated limit, 

lim lim \ = T. 

By means of the numerical blow-up rate we observe a propagation property for 
blow-up points. We prove that the nodes adjacent to those that blow-up as the 
maximum may also blow-up (opposite to the continuous problem), but they did it 
with a slower rate and the number of adjacent blow-up nodes is determined only 
by p and is independent of h and A. 

In other words, if we call B*{U) the set of nodes k such that 

lim ui(T„,A-t^y/(^-i) - Cp, 

j-»oo 

the number of blow-up points outside B*{U) depends (explicitly) only on the power 
P- 

We split the paper in two parts, in the first part we deal with the explicit scheme 
and in the second one we develop the analysis for the implicit method. 



2. The explicit scheme 

The main tool in our proofs is a comparison argument, so first of all we prove a 
lemma which states that this comparison argument holds. Since we need restrictions 
in the time-step to prove this lemma they are essential for every result in this section. 
That is not the case of the implicit scheme. 

Definition 2.1. We say that {Z^) is a supersolution (resp.: subsolution) for Q if 
verifies the equation with an inequality > (<) instead of an equality. 
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Lemma 2.1. Assume the time step verifies 

r," < mm — . 

i<i<N an 

Let (?/"'), {IP) o- super and a subsolution respectively for (0)) such that < if , 
then IP < [/"' for every j . 

Proof: Let = [/"' — IP, by an approximation argument we can assume that we 
have strict inequahties in Q, then {Z^) verifies 

/-.x ii/az^+i > -Az^ + M{{W)p ~ {iP)p, 

^ ' Z" > 0. 

If the statement of the Lemma is false, then there exists a first time and a 
node Xj such that z'i'^^ < 0. At that time we have 



z 



> (1 - r,!^)zi + r, - ^ a.,zi + -&\>Q, 
a contradiction. □ 



2.1. When does the solution blow up. In this section we find conditions under 
which the solution of (O blows up, we begin with some lemmas. 

As the matrix A is a symmetric (property (P2)), there exists a basis of eigenvec- 
tors for the following eigenvalue problem 

A(j)i = \M(j)i. 

We call rj = ri{h) the greatest eigenvalue of this problem, that is 

< Ai < T]{h). 
Lemma 2.2. For every y S M.^ there holds 

{Ay,y) <7j{h){My,y). 

Proof: As the matrix M defines a scalar product in M.^ , we can assume that the 
eigenvectors (jfi are normalized such that 

Let y e R^, y = J2iLi o^i4>i^ then 

/ N N \ 

{Ay,y) = (^a,X,M(f>,,^aj(t>j\ 

\i=l 3=1 I 

N 

1=1 

< r]{h){My,y). 

□ 

Lemma 2.3. If {U^)j>o is large enough, then blows up infinite time. 
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Proof: Recall the definition of 



N 



= ^ ruiui , 

i=l 

and observe that there exists constants c, C > that depend only on h such that 

c max u] <w^ <C max u]. 

l<k<N ' l<i<N ' 

Now we observe that we can choose j > jo in order to get large enough to 
verify 



N N N 



k=l i=l i=l 



=1 i=i 

AT 



2 

i=l 



> + cTj {w^ y 

= + cA. 



Applying this inequality inductively we obtain 

> + cX{j — jo) > cj, 

hence — > oo as j ^ oo, it remains to check that Xi'^i < ^o do that we 
observe 

A ^ A 



^ {wi)P ~ (wO + jcA)f ' 
and 

^ A t°° ^ . 

y -r-F, < / T-^ as < 00. 

fr{iw°+j (cA))P - Jo {w° + csX)P 

This completes the proof. □ 

Remark 2.1. In the course of the proof of the above Lemma we showed not just 
that blows up, we also proved > cj. 

Now we are going to prove the reverse inequality to obtain the asymptotic be- 
havior of ||?7-' ||oo- 

Lemma 2.4. If (U^) is unbounded then 

WU^Woo-^W^ ^j 

Proof: The relation \\U^\\ ~ w-' is trivial since they define equivalent norms in 
R^. So we just have to prove < Cj. To do that we observe that 
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N N 



N 



J + 1 



< W-' 



k=l 1=1 
N 



i=l 

= + ex. 



We apply this inequality inductively again to get 

< w" + CXj < Cj, 

as we wanted to prove. □ 

Theorem 2.1. Assume the time step tj verifies the restriction 

2 

Then positive solutions of (0 blow up if there exists jo such that ^fi{U^'') < 0. 

We remark that the condition ^ii{U^°) < is similar to the one for the blow-up 
phenomena in the continuous problem, in fact for the continuous problem this is a 
necessary condition. 

Proof: First we observe that <i>/i decreases with j, in order to do that we take 
inner product of ^ with C/^+^ — to obtain 

= {—M{W+^ - W) + AW - M{WY, - W) 

+ {AW, - ^{AU^, U') - i(Mp(^Jf -\ (t/^+i - U^f). 
Hence we obtain, 

< Tj-(r /^"''^'^'' ' - \){MdU^+^,dU'+^) + ^{AdW+\dW+^) 

due the restriction in the time step tj and Lemma lT^ Actually < ^fi{U^) 

unless {U^) is independent of j. So, is a Lyapunov functional for Q). Next we 
observe that the steady states of (@J have positive energy. Let (W^) = be a 
stationary solution of then 

= ~AW + MWP. 
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Multiplying by we obtain, 

= -hAW,W) + P^^ [ lMWP,W) 
2 2 p + 1 

> -<^h{W). 

Now, assume (U^) is a bounded solution of then there exists a convergent 
subsequence of (U^) that we still denote (U^). Its limit is a steady state with 
positive energy. 

As ^h,\{U^) decreases and there exists jo such that ^h.xiU^") < then ^h,\{W) < 
0, a contradiction. We conclude that ([/■') is unbounded and by Lemma [2.31 has 
finite time blow-up. □ 

Corollary 2.1. Assume the time-step restriction of the above theorem and the 
convergence hypotheses (HI), (H2). Let uq an initial data for such that u blows 
up in finite time T. Then u^^x blows up in finite time Th \ for every h, \ = X{h) 
small enough. Moreover 



lim lim x — T 



Remark 2.2. If the fully-discrete method converges in Hq{Q?) a.e. t then A can be 
chosen independent of h. 

Proof: We observe that if u blows up in finite time T then (see |15| ) 



Jn 2 Jn 

This implies that there exists a time t^ <T with $(u)(to) < 0. Now we use the 
convergence of u^^x to Uh in [0,to] and the convergence hypothesis (HI) to see that 

lim lim ^h.\{uh.\)[ta) = ^{u){to). 

h^O X^O 

So for h,\ — X{h) small enough we get ^h,x{uh,x){to) < and so, by the above 
Theorem {U-') blows up. 

Now we turn our attention to the blow-up times. In it is proved that the 
blow-up time of the semi-discrete solutions (solutions of ) , that we are going to 
denote Th, converges, as /i ^ 0, to T. In that work the poof is given for Q = (0, 1) 
and a finite element method. Despite for this case the proof is very similar, we 
sketch it for the sake of completeness. In jj^ is proved that if the continuous 
solution blows up then for every h small enough the semidiscrete scheme also does. 
Hence we can assume that the semidiscrete solution U{t) is large enough in order 
to verify 



j^{MU{t),U{t))^2{MU'{t),U{t)) = 
2{-AU{t), U{t)) + 2{MUP{t), U{t)) = 
-A^h{U{t)) + ll(Af(7f (t), U{t)) > 
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Integrating between to and Th we obtain 

(3) in - to) < — —■ 

where C depends only on p. 

Given e > 0, we can choose M large enough to ensure that 

C \ e 
< -. 



2' 



As u blows up at time T we can choose t < | such that 

-$(u(-,r-r) > 2M. 

If h is small enough, 

-^h{U{T~T))>M, 

and hence 

n-{T-r)< ( ^ —] < (^) < 

Therefore, 

\Th-T\ < \Th - (T - t)\ + \t\ <e. 
We have proved lim^^o Th = T, so we just have to prove that for fixed h 

lim Th,\ = Th. 

A-+0 

To do that we observe that from Lemma 12.31 we know that there exists jo, that does 
not depend on A such that for j > jo 

hence 

A 



^ A ^ p A 

" + cA(Z - jo)Y - J J [wio + cA(s - jo))P 

= — / — ds < — — as. 

C Ju]io+c\(j-jo) ^ Jwio SP 

(4) 

This holds for any jo large enough and for every j > jo . In particular we get 

1 r°° 1 

Th.x -t' <- —ds. 

c J^i sP 
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This inequality has a great meaning. It says that if is large, then P is close 
to blow up (independent of A). So now, given e > we can choose K large enough 
in order to have 

- —ds<-, K P < -. 
cJk sP - 3 

Next we choose r < | such that 'Y^mkUk{Th — 2r) > 2K (remember that 
(Mi(i), . . . , Uff{t)) is the solution of the semidiscrete scheme). For A = X{h, t) small 
enough we get, from (H2), that > K for every j such that — 2r < P < T^ — t. 
We choose one of those j and compute 



|r/,,A-T,,| < \Th,x^P\ + \P -Th\ 

< - / —ds + 2T 

< £ 

□ 

2.2. Blow-up rate. In this section we study the asymptotic behavior of numerical 
solutions with blow-up. 

Theorem 2.2. Let Uh,\ a solution with blow-up at time Th,\, then 

max u] - {Thx-t'y^'^P-^\ 

l<i<N ' ' ' 



Moreover 

,i/(p-i) 

lim max 

j^oa l<i<N 



We want to remark that this is the behavior of the continuous solutions with 
blow-up. 

Proof: We know from Lemma l2 . 31 that — ^ miU'f verifies 
so we have 



k=j + l k=j+l ^ ' 

~ J, {^{s))P 

(5) 

Here w{s) is the linear interpolant of [w^ = w(j)), hence for j < s < j -I- 1 we have 
w'{s) = w-'"'"^ — > cA, and so 
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or equivalently 



max ui < Cw^ < C(Th x - t^)-^/^P-^'> . 

l<i<N ' V , / 



The inverse inequalities can be handled in a similar way to obtain 



max u^' ~ ~{Thx- t^Y^'^^-^^ 

l<i<N ' \ , / 



Next we recover the constant Cp in the asymptotic behavior of the numerical 
solution, to do that we will change variables but first we prove a short lemma. 



Lemma 2.5. 



Proof: 



lim^i^ = l 



□ 

Now we change variables, in a way inspired by ^Jj^Hl- Let {Y^) be defined by 

(6) yf =ui(T^^^^t^y/(P-^) l<z<N. 

In the sequel of the proof we will use ^yf^^ to denote 

i+i i 

yj -yj 

T,i{Tu,x-py 

This can be thought as Tj/{Th^x — t^) to be the time step in the new variables. 
With this notation the new variables verify 



{Th,x-P)— tt 



(7) {n,x-p)- 



[in.x-t^+')^ -{Th,x-tn^) 



y° = iTh,xy/^P-^'>uoix,), l<i<N+l. 

Now assume there exists e > and a subsequence that we still denote {yf ) such 
that yf > Cp + e for some i — Then for those yf we have 
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We also know from the blow-up rate that the new variables are bounded and 
so, applying Lemma 12.51 we obtain for j large enough 



(8) 
(9) 

(10) 



> 



p-V 



- 1 



> 5. 



This means that actually yl > Cp + e for every j large enough and consequently 
(O is verified for all those j. So yf is unbounded, a contradiction. 

If we assume < Cp — e arguing along the same lines we obtain that yj verifies 



miAyf~^^ < S + m. 



ivi) 



p-V 



1 



{Th,x-v+^)— {n.x-^^)—' 



< 26 + C[{ylY 



p-l 



(11) 



This shows that either yl — > as j ~> cxd or niiAyj < —S, which means that 
yj is not bounded from below (this is not possible). 

We conclude that if yj does not go to zero, then it goes to Cp, as the blow-up 
rate implies that for every j 



we have 



max yl > c, 

l<i<N ' 



lim max yl = lim max (Th x ~ t-'Y^^^ ^-'u^ = Cr,, 

j~^ool<i<N * j^ool<i<Af' ' ' IF 



as we wanted to prove. 



□ 



2.3. Blovir-up set. Now we turn our attention to the blow-up set. In order to 
do that we consider the set B*{U) composed of those nodes that blow-up as the 
maximum (a precise definition of B*(U) is given in the introduction) and we study 
the behavior of the adjacent nodes, then we repeat the procedure with these nodes. 
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Definition 2.2. We define the graph with vertices in the nodes and say that two 
different nodes are connected if and only if aij ^ 0. We consider the usual distance 
between nodes measured as a graph, see [TBj . Finally, we denote by d(k) the distance 
of the node Xk to B*{U) also measured as a graph. 

We prove that blows up if and only if d{k) < K where K depends only on p, 
Tiieorem 2.3. Let B*{U) be the set of nodes, {xk}, such that 

ui^{n,x^t^y^. 

Then the blow-up propagates in the following way, let p > 1 and if G Nq such that 

< P < -^^^ (K is the integer part of l/{p — 1)). Then the solution of 
blows up exactly at K nodes near B*{U). More precisely, 

ul +00 d{k) < K. 

Moreover, if d{k) < K, the asymptotic behaviour of (w],)j>i is given by 

ui^{n,x-t^r^+''^''\ 

^fP + ^ and ifp = d{k) = K 

ui^\n{Th,x-t^). 

We want to remark that more than one node can go to infinity, but the asymp- 
totic behavior imposes ^ — > (j — ^ oo) if d(k) > d{i). 

Proof of Tiieorem 1.4 We want to show that the blow-up propagates K nodes 
around B*{U), we begin with a node Xk such that d{k) — 1. We claim that the 
behavior of uj, is given by 

( j-P+^ if p<2 
~ < Inj if p = 2 
[ C if p > 2, 

to prove that we will show that 

J 

^ A^STs-1, 

3 = 1 

which has the behavior described above, can be used as super and subsolution for 
an equation verified by choosing A appropriately. 

We observe that ul. satisfies 

N 

mkdul^^ = - ^ aikuj + mk[vPf.Y 
1=1 

^ Ci( max ul) - C2ui + Cs,mk{u'uY ■ 

l<i<N " " 

This means that there exists constants Ci > 0, i = 1, 2, 3 such that for j large 
enough 

(12) < c,j - cWk + C3(ui)f 
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and 

(13) dui+'>c,j-C2ui + c,{uir 

Now we observe that if A and j are large enough, verifies 

dw^^ = A{j + 1) 

since [w^Y / j Q as j goes to infintiy. Hence is a supersolution for l(T^ and 
so 

On the other hand if we choose A small we get 

dw\ = A{j + 1) 

< Cij - C2W^ + C3(w^)^', 

Hence now we can use as a subsolution for (|15|) to handle the other inequality. 
Therefore 

ui - w\. 

We observe that if p < 2 the node Xk is a blow-up node and we also have the 
blow-up rate for this node (u^ ~ J"^"*"^)- If p > 2 this node is bounded. Next 
we assume p < 2 (if p > 2 it is easy to prove that every node k with d{k) > 1 is 
bounded) and we are going to find the behavior of a node, that we still denote fc, 
such that d{k) = 2. That is, it is not adyaccnt to B*{U) and it is adyacent to a node 
wich has the behavior j^^'^^. 



s=l 



Now let 



and observe that verifies 



N 
i=l 

This means that there exists constants Ci,Ci > 0, z = 1,2,3 such that for j large 
enough 

(14) < - C^2< + C3(«i)P 
and 

(15) dui+^ > - C2ui + C3{ui)P 
Now we observe that if A and j are large enough, w\ verifies 

dw\ = A{j + l)-f+2 
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since {w\Y / j as j goes to infintiy. Hence w\ is a supersolution for l(T^ 

and so 

On the other hand if we choose A smah we get 

Now we can use w\ as a subsolution for l(T5|l to handle the other inequahty. So 

We observe that if p < 3/2 the node Xk is a bfow-up node and we also have the 
blow-up rate for this node (it^ ~ j '^^'^^)- If p > 3/2 this node is bounded. In the 
case p < 3/2 we repeat this procedure inductively to obtain the theorem. □ 

3. The implicit scheme 

In order to avoid the time step restrictions we now introduce an implicit scheme 
and prove that similar properties can be observed. However we have to remark 
that near the blow-up time the adaptive procedure forces the time step to be much 
smaller than the space discretization parameter h. This suggest that an adequate 
method could be to begin at time zero with the implicit scheme in order to avoid 
time-step restrictions and, as when the solution increases the time step is reduced, 
if the solution has finite time blow-up, then there will be a moment such that the 
time step restriction for the explicit scheme will be verified. From that time one 
can continue either with the explicit or with the implicit scheme. As in the explicit 
scheme section we begin with the comparison Lemma. 

Liemma 6.1. Let {U'), (LP) a super and a subsolution respectively for ® such that 
< , then IP < for every j . 

Proof: Let ^U-W, we assume that we have strict inequalities in ©, then 
(Z^) verifies 



f.. MdZ^+^ > -AZ^+^ + Mi{U^)P~{lP)P), 

^ ' ZO > 0. 

If the statement of the Lemma is false, then there exists a first time and a 
node Xj such that zf"*"^ < 0. There we have 



-^+' - - ->^4+^ - —4^^ + (^r - (Mir > 0, 



a contradiction. 
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3.1. When does the numerical solution blow up. 

Lemma 3.2. If {U^)j>() is large enough, then blows up infinite time. Moreover 

\\W Woo - r^j 



Proof: 



N N N 

-LI 

fe=l i=l 
N 



< + Tj mi{u{Y 

i=l 

< W^+CTj{w^)P 

= + ex. 



Hence < Cj. To prove the inverse inequality we observe that 



N N N 



k=l 1=1 i=l 

(2) 

that is 



(3) (1 + ClTj>J'+l > + C2Tj{w^)P. 

Now we look for a subsolution of Q of the form ~ Tj. This sequence verifies 

(1 + Cirj>J+i = + rCiTjj + r(l + Cir,) < z^ + C2T,iz^)P 

if r is small enough. As the discrete maximum principle holds for this equation we 
obtain 

> Z-' = Tj. 

This completes the proof. □ 

Theorem 3.1. Positive solutions of |(SJ) blow up if there exists jo such that<^h{U^°) < 
0. 



Proof: Again we first observe that ^h{U^) decreases with j, to do that we take 
inner product of © with J7^+^ — to obtain 
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Hence we obtain, 



< 0. 

This implies that $h is a Lyapunov functional for ©. We observe that the steady 
states of (jni are the same of Q, so they have positive energy. Now, assume (U^) is 
a bounded solution of then there exists a convergent subsequence of ([/■') that 
we still denote (U^). Its limit 14^ is a steady state with positive energy. 

As ^h{U^) decreases and there exists jo such that ^h{U^°) < then ^h,\{W) < 
0, a contradiction. We conclude that {U^) is unbounded and by Lemma [3.21 has 
finite time blow-up. □ 

Corollary 3.1. Assume the convergence hypotheses (HI), (H2). Let uq an initial 
data for such that u blows up in finite time T. Then Uh.x also blows up in finite 
time Tfi^x for every h, \ = \{h) small enough. Moreover 

lim lim Th \ = T. 
Proof: If u blows up in finite time T then (see Ill l -llS p 

Jo 2 P + 1 

Hence there exists a time to < T with <!>(«) (to) < 0. Now we use the convergence 
hypothesis (HI) and the convergence of Uh,\ to Uh in [0, to] to see that 

lim lim ^h.x{uh.\){to) = ^u){to). 

So for h, X{h) small enough we get ^h.x{uh.x){to) < and so (U^) blows up. The 
convergence of the blow-up times is obtained like in the explicit scheme. 

□ 

Next we turn our attention to the blow-up rate of the discrete solutions. 
3.2. Blov^f-up rate. 

Theorem 3.2. Let Uh,x a solution with blow-up at time Th,x, then 

max ui^in.x-t')-'^^"-'^. 

l<i<N 
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Moreover 

/ 1 \i/(p-i) 
lim max ui{n,x - t^y/(P~^^ = Cp = { 

i^oo l<i<N * ' ' \p — 1 J 

Proof: The first part of the proof is the same as the one for the explicit scheme so 
we assume we have proved 

and we are going to recover the constant Cp. We change variables as in the explicit 
scheme, let (Y^) be defined by 

(4) yj^u>(Th,x-t^y/^-''^ l<i<N. 

In the sequel of the proof we will use Ayf'^^ to denote 

J + l 3 

yj - VI 

this can be thought as Tj/{Th,\) to be the time step in the new variables. With 
this notation the new variables verify 

mi^yi = -{Th,\-t^)2_^ak,yi +^i-Ff _p\i/(p-i) ^yil 

1=1 ^ ^''^ ' 

(5) 
(6) 

If we assume that there exists e > and a subsequence such that yj > Cp + e 
for some i — Then for those yf, as they are bounded, if j is large enough we 
have 

Ayl+' > 



> 

(7) 

This means that actually y^ > Cp + e for every j large enough and consequently 
(jZj) is verified for all those j. So yj is unbounded, a contradiction. 



p-l 



{Th X - t^'+i)i/(P-i) 



p-V 



5. 
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The case where there exists an infinite number of yj such that yj < Cp — e is 
very similar. So we can conclude that as ^ cxd either yj — > or ~* Cp. Now 
we use the blow-up rate to obtain 

lim max = lim max (Th x - t^V^'-^^'^^ui = C„, 
as we wanted to prove. □ 

3.3. Blow-up set. The propagation property for the blow-up nodes holds for the 
implicit scheme and its proof is very similar. We do not include it. 



4. Appendix 



In this appendix we prove that if the general method considered for the space 
discretization is consistent (see below) then the totally discrete method converges in 
the L°° norm. We perform the proofs for the explicit scheme, they can be extended 
to the implicit one. 

Definition 4.1. Let w be a regular solution of 

wt^ Aw + f{x,t) mnx(0,T), 
w^O ondVtx{Q,T). 

We say that the scheme ((5J is consistent if for any t £ (0, T — t) it holds 

N 

(1) niiWtixi.t) = - ^ aikw{xk,t) + mif{xi, t) + pi^t), 

fe=i 

and there exists a function p : — » such that 

max ^^''''^ -^^ < pih)^ for every t G {0,T — t), 
i rrii 

with p[h) — > if h ^ 0. The function p is called the modulus of consistency of the 
method. 



If we consider for example a finite differences scheme in a cube Q. = (0, 1)'' C M'^. 
Then the modulus of consistency can be taken as p{h) = Ch^ . 

Theorem 4.1. Let u he a regular solution of ^ G C^'^(ri x [0,T — rjj anc? u^.a 
the numerical approximation given by (0 then there exists a constant C depending 
on \\u\\ in C^-^{n x [0,T-r]) such that 

\\u - M/i,A||L~([0,T-r],L2(O)) < C{p(h) + A). 

Proof: 

We define the error functions 

4 = u{xi,tj)-ui. 

By these functions verify 

N 

mj9e^^^ = - ^ Oikel + mi(u''{xi,tj) - {u{Y) + pi{h) + CniiX, 
fc=i 
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where C is a bound for ||utt||L«=(Ox[o.T-T])- Let 

to — max{i : t < T — t, maxmax |e^| < 1}. 

i tj<t 

We will see by the end of the proof that = T — r for h, X small enough. In [0, ip] 
we have 

N 
k=l 

hence, in [0,<o], = (ej,...,e^) satisfies 

(2) MdE^+^ < -AE^ + KME^ + {p{h) + CX)Ml\ 

Let us now define — (w^, . . . , WAr(i)), which will be used as a supersolution. 

wl^e^^^+^^'^{\\em^ + p{h) + CX). 

It is easy to check that verifies 

MdW'+^ > -AW^ + KMW^ + {p{h) + CX)Ml\ 

here K is the Lipchitz constant for f{u) — in [0, T — t)\\i^-x]. Hence is 
a supersolution for and by Lemma [2. II we get 

4 < e('^'+^^*^(||e°|Uo.(o) +PW + CA), e [0,to]. 

Arguing along the same lines with —E^, we obtain 

H\ < e(^^+i)^(||£;(0)|U + p{h) + CX) < C{p{h) + X), t, e [0,to], 

by our hypotheses on the convergence of the initial data. Using this fact, since 
p{h) goes to zero, we get that |e^ | < 1 for every tj G [0, T — r] for every h, X small 
enough. Therefore = T — t for /i, A small enough. This proves the convergence 
of the scheme. In fact we have that for every h < ho, X < Xq 

max max — u{xi,tj)\ < C{p{h) + A). 

j l<i<N 

□ 
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